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By focusing on the interchangeable role in a generating function (i.e., j3 <-> E in the Laplace 
transform), the superstatistics proposed by Beck and Cohen can be viewed as a counterpart of 
the canonical partition function. Some formal properties of this superstatistics are presented in 
connection with thermodynamic structures and information aspects. For any combination of the 
local equilibrium statistical factor and the form of fluctuating field, which are ingredients of making 
a generic superstatistics, a variance of the fluctuating quantity appears in the correction term. This 
fact enables us to relate parameters contained in the statistical factor with the variance in any 
situation. 

PACS numbers: 05.90.+m, 89.90.+n, 89.20.-a 
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I. INTRODUCTION 



Superposition of statistical factors can be considered as 
a new attempt to try to formulate a statistical mechan- 
ical description for unusual systems such as ones having 
a fluctuating temperature field. Since, in the traditional 
equilibrium statistical mechanics, we suppose that sys- 
tems are put in constant with a heat reservoir (in canon- 
ical ensemble) , situations that arise inhomogeneity in in- 
tensive variables over spaces considered may be outside 
the scheme of conventional treatments^, 

In this paper, we will be particularly concerned with 
presenting some detailed structure of the notion of the 
superposition of statistical factors and its relation to the 
fluctuation. The purpose of this paper is to understand 
what aspects of a system the superposition of the statis- 
tical factors delineates. To this end we especially take 
up of a mixture of the Boltzmann factors, which is re- 
cently proposed by Beck and CohenQ- They named this 
quantity superstatistics. All considerations through this 
paper is a situation where local inverse temperatures (3, 
which is usually assumed to be a constant in the ordinary 
canonical ensemble theory, are not homogeneous in space 
but differ by positions. 

The superstatistics denoted by symbol B(E) allows the 
infinite types of system's distribution in energy E, once 
the fluctuating distribution f{(3) is given. Three crucial 
premises taken by Beck and Cohen are (a) a system is 
partitioned into cells that can be considered to be reached 
an equilibrium locally, which is characterized by a single 
intensive variable (b) its statistical factor is Gibbsian (c) 
the separation between two time scales is adequate, that 
is, the time for relaxing to the each local equilibrium 
state is much faster than that for varying the distribu- 
tion of (3. Under these assumptions B{E) is calculated 
by integrating f(f3)e~ l3E with respect to (3 from zero to 
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infinity. Then the usual Boltzmann factor for the system 
is obtained only when the system has homogeneous tem- 
perature /? = (J3) i.e., f(f3) = 8(j3 — (/?)). In this view, 
the linear average of the Boltzmann factors over (3 can 
contain information that characterizes the macroscopic 
behavior. 

From a formal point of view, Beck and Cohen ex- 
emplified that the correction to the equilibrium factor 
e -(P)E universally follows (1 + a 2 E 2 /2), as long as oE is 
small, where a is the standard deviation of f((3). Espe- 
cially when f((3) represents the y 2 distribution the asso- 
ciated superstatistics produces the Tsallis type power-law 
distribution 0, Q, and by comparing with the correction 
term (1 + (q - l){(3) 2 E 2 /2) for the x 2 distribution, the 
parameter q contained in the statistical factor can be 
formally related as q = 1 + (o-/((3)) 2 . Then the case q = 1 
corresponds to no fluctuation in /?, i.e., constant temper- 
ature over the system. Once the B(E) is given, the forms 
of entropy and the associated energy constraints can be 
determined[6( and the maximum-entropy principle can 
be a good tool for estimating a statistical properties 

As of today, we have taxonomy for the applications of 
superstatistics (inexhaustive list): the log-normal type in 
turbulence M, El and the x 2 type in solar flares 01) cos ~ 
mic ravs|ll|. and wind velocities in airport 12]. A violent 
gravitational process can be another possible class of su- 
perstatistics, where the Gaussian distributions in velocity 
of stars are equally superposed by temperature Il4j . 
Other recent applications of the superstatistics are briefly 
reviewed in [lfij . 

At the present status, however, there seems no guiding 
principle (dynamical mechanism) to choose any specific 
combination of the forms for the local distribution and 
the fluctuation distributions. We must only refer to the 
mechanism as a certain complex dynamics and can not 
answer how complex it is. In this sense, the present sta- 
tus on its necessity for introducing B(E) may still be 
contentious^^, 0|. This ambiguity is also related to the 
fact that some phenomenological data can be fitted for 
reasonable regions with a generalized exponential func- 
tions having more than two parameters |l£|, and this fact, 
however, does not imply that the system can be uniquely 
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described by the underlying generalized statistical me- 
chanics based on the parameterized entropy which should 
be maximized. 

These facts are our main reason for avoiding going into 
discussion of a specific physical system in this paper, in- 
stead we shall explore some formal connections to the 
thermodynamic relation and to information structure. 
The plan of this paper is as follows. In Sec. II, we in- 
vestigate six properties possessed by B(E), which have 
not been presented so far. The subsections A and C are 
concerned with the connection to thermodynamic struc- 
tures. From a view point of information, we provide a 
universal relation by considering one dimensional system 
in subsections D. In subsections B, E and F we give re- 
sults which are direct consequences from the definition 
but are worth presenting in future developments. Sec. Ill 
devotes to a discussion for the notion of superposition of 
statistical factors and to a summary. 



which becomes F{(3) — > oo(/3 — > oo), we write 



V>(£) 



- (3E dF(f3). 



(4) 



Then we can evaluate this as 



1>{E) > e bE / e^ E dF(/3) > ^-E^p+b)^ 

J min/3 

_ e E(b-rninf3)/2j /g) 

where I = f£^ff +l ' dF((3) is a finite. Therefore, we 
have i/)(E) — > oo when E — > oo. Since the second deriva- 
tive of \ntp(E) with respect to E is the same as that of 
\nB{E), the property that tp(E) -> oo(£; -> 0,E -> oo) 
together with the property (iii) concludes that there 
should be a single minimum value E* at which 



II. SOME FORMAL PROPERTIES OF B(E) 
A. Existence of average in /3 



dhxip(E) 



dE 



E=E* 



, B'(E*) n 



is realized. 



(6) 
□ 



Let us assume that the system considered is situated 
in an environment where few partitioned region has large 
temperature and the number of regions in the system 
grows steadily as their temperatures go down. The most 
probable distribution for the regions would be the in- 
creasing function of (3. From the definition, the super- 
statistics has the following properties: 

(i) B(E) -> oo when E ->• 0. 

(ii) B(E) is a C°°-class function with respect to E for 
any E > 0, i.e., 



B^(E) = (-1)" / dl3f3 n f(f3)e-? E , n = 1, 2, . 

where we set the minimum in (3. 

(iii) \xiB(E) is a convex function for E > 0. 

This fact follows from the direct calculation below: 



(1) 



dE 2 



In B(E) = 



B{E)B"{E)-(B'(E)f 



B 2 {E) 



B{E) j 

min(3 

> o. 



(2) 



As an important consequence of the above properties ((i)- 
(iii)), we can prove the following statement: 
For b which satisfies < min(3 < b, there exists one single 
positive solution for any E > such that 



-dE lnB ^ 



(3) 



Proof. Putting ip(E) — e hE B(E) and using the integrated 
density of states F(j3) for the temperature distribution, 



Eq.Q can be regarded as a counterpart of the fun- 
damental relation — din Z/d/3 — E, which ensures the 
existence of energy in the system in usual statistical 
mechanics We shall find that b represents the av- 
erage inverse temperature in section Til El 



B. Behavior with low temperature cut-off 

Next, we are concerned with a behavior of the super- 
statistics at higher temperature. If /(/3) has an asymp- 
totic approximation, 

f(J3) ~ a /? 70 + ai/? 70 + • • • + a„/3 7 " , as -► oo (7) 

where exponents are assumed to to be satisfied with — 1 < 
7o < 7i < • ■ ■ < 7n- The superposition of the factor with 
the low temperature cut-off max[3 is defined 



B cut (E) 



maxf3 



f(f3)e-P E df3. 



(8) 



Note that B cut (E) — > B(E) when max {3 — > oo. Then the 
Watson's lemma 201 states that 



B cut {E)~Y^CLkT{lk + l)E-^\&s E^w. (9) 

fc=0 

This states that when f{(3) consists of some different 
power laws in (3, the corresponding superstatistics is ex- 
pressed by the superposition of the power laws in energy 
with different exponents. 
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C. Representation by specific heat 

In this section we delineate the superstatistics in the 
light of a thermodynamic quantity i.e., the specific heat 
of the system. As was supposed in section III Al we as- 
sume that f(0) represents the density of states of inverse 
temperatures in a system around 0, which monotonically 
increases with 0. Then the probability to find a system 
within the interval and (3 + d/3 can be expressed as 



p(0)d0 oc e- 0E f(l3)d(3. 



Then the integrand e lnB +0 E tells us that a quantity 
4>(0,E) = h\B(E) + 0E provides information on sys- 
tem's fluctuating field. Note that this is the same quan- 
tity by the Legendre transform of B(E), which will 
be introduced later in this section. Also, the reader 
should refer to a discussion revolving around the Laplace 
methods pll] for extreme cases E — > oo,0 — » and 
E -» 0.3 -» oo[l7l 15^. We instead focus on the 
evaluation for the intermediate energy region. At a 
saddle point (E ) of the energy coordinate, <j>(/3,E) 



(10) should satisfy <j>' = (d<f>(J3,E)/dE) E=Ea = and 



Since the Boltzmann factor decreases with 0, p(0) should 
have an extremum at a certain value of — 3*, so that 



(d z 4>(3,E)ldE 2 ) E =E > 0. Then when we expand 
4>{0, E) around E = Eq, we have the expression 



d_ 

80 



[e-^f(0)} 







0=0* 



that is 



08 



[In /(/?)] 



= E. 



(11) 



(12) 



<t>{3,E) ~</> + -(E - E Q ) 2 <b' ' 
where </> — <fi(8, Eq). Therefore we can write 



0=0* 



Expanding \np(0) around 0* 



ln[p(0)} = (-/?*£ + In f(0*)) 



i a 2 



-0E 



2d0 2 
-0*E + ]nf(0 



0=0* 



2k B 



(0-0*? + 
(0-0*f. 



(13) 



where C v — (dE/dT) is the specific heat. Thus we can 
evaluate B(E) as 



B(E) = / P(0)d0~e 



-0*E+\nf(0*) 



yf2Trk B 0*/C v e-l 3 * E + hi M'\ 



d,0 
(14) 



Since the specific heat can be usually considered as be- 
having O(N), N being the number of particles in the sys- 
tem, the term which contains C v contributing as 0(ln N) 
in hi B(E) becomes negligible in comparison with others 
for large N. In that case —0*E + In /(/?*) is mainly re- 
sponsible for the thermodynamic property of the system. 
However we must take into account C v for systems with 
small N. 



D. B(E) as an information and its correction by 
fluctuation 

To illustrate more information aspects on the super- 
statistics by fluctuation f(0), let us take a look at another 
point. The fluctuation distribution can be calculated by 
an inverse Laplace transform with keeping E constant: 



-1 />C+200 

f(P) = 77- / B(E)e? E dE. 



(15) 



2ttI 



c+«oo 



e (E-Eo) 2 <P /2 dE ^ 



(16) 



(17) 



A change of variables as c = Eq and ix = E — Eq (x e R) 
provides an expression 



\/ 2lT( t>'a' 



(18) 



Thus, the corrected bit 7wm&er|23j of the system at equi- 
librium, which is defined as ln/(/3) = Imc{0), incor- 
porated small fluctuations around equilibrium energy E, 
reduces to 



lmc 



09) 



1 



In <j> 



(19) 



This plays a main role in this study and reminds us the 
application to the canonical black hole entropy by Das et 
aL[24|, where the general corrections to the microcanoni- 
cal entropy of some black holes were calculated. The suf- 
fix MC is intended to denote macroscopic because f(0) 
has a global information for the system. (j> (0, Eq) in the 
above equation is found to be equivalent to (0 2 ) s — {0)1 
if we employ the average of and 2 by means of p s as 
we will see in the next section, which allows us to express 
as I MC (0) ~ O - \n[(0 2 ) s - (0) 2 s }/2. 

In order to delineate further structure which is uni- 
versally possessed by B(E), we consider one dimensional 
system whose segment at position n is assumed to be 
in a local equilibrium prescribed by a Boltzmann fac- 
tor e~ /3nS . In order to see the effect of the fluctuation 
of temperature on the bit number of whole system, we 
can borrow an approach that was taken by |25| for the 
canonical entropy of very large black holes obeying the 
Beckenstein-Hawking area law|2^. In this approach, it 
was shown to have the universal logarithmic correction to 
the microcanonical black hole entropy. The superstatis- 
tics for the present system may be expressed as 



B(E)= J2 f(^)e 



-PnB 



(20) 



4 



We would be appropriate to interpret this form as 
weighted sum of the local factors according to the value 
of inverse temperature of the segment n. We can move 
to a continuous expression 

/oo 
f(l3{x))e-^ E dx, (21) 
-CO 

where we have applied the Poisson resummation formula, 



E 



E 



m— — oo 



dxg(x)e 



(22) 



which holds for sufficiently large x with ^ TO e ~ l27Tmx ~ 1. 
Then by converting to the integration with respect to /3 
we obtain an expression 



B(E) 



exp[I MC (f3) ~(3E- In (23) 



By employing the method of steepest descent for the sad- 
dle point /3q, we have 



B(E) 



2tt 



x exp 



hic(M-PoE- 



ln 



(d/3\ 




\dx) 


0=00, 



(24) 



Here we define Ic (P) by the Legendre transform of pairs 
(E,-\nB(E))^((3,IcW),le., 



I c {/3)=lnB(E)+(3E, 



(25) 



which is the same as 4>{f3, E) introduced in the begin- 
ning of this section. This is an analogy of the standard 
thermodynamic definition of canonical entropy S c from 
the canonical partition function Z c \ — \nZ c = —S C + (3E. 
Then Ic(fi) can be interpreted as a quantity which con- 
tains information of fluctuation and spacial inhomogene- 
ity. Therefore Ic(P) evaluated at the saddle point (3 — (3q 
is calculated to be 



Iciflo) = WA,)-ln(j| 

1, ,7T N 

+ 2 ln ^- 



by 



0=00 



(26) 



AI = — In 
2 



7T dp 2 



Imc(Po) 


= AI 






\dx ) 


0=00 



(27) 



This formula is expected to hold for the most general 
circumstances without depending on particular choice of 



/(/3) and (3{x). However observe, for example, that with 
a constant temperature gradient in x, which is a possible 
setting for systems, the quantity in the square brackets in 
Ea. (|27|) can be negative for some /(/3)[27|- Moreover, a 
uniform distribution and finite level distributions always 
give d 2 lMc/d(3 2 = 0. This explains that the saddle-point 
approximation for the calculation of the effect of fluctu- 
ation (the logarithmic shift term) to the macroscopic bit 
number is limited to the case where I MC ([3o) < is sat- 
isfied. 



E. Discrete levels for B 1 (E) without quantization 

In this section, we seek yet another connection to ther- 
modynamic expression and the superstatistics. What we 
are concerned is a relation with the variance of (3. To see 
this, we introduce a conjugate distribution function |l9j| 



-/3E 



B{E) 



p s d(3 



(28) 



Then 



d 2 \nB{E) _ B"(E) ( B'(E) Y 
dE 2 ~ B(E) \ B(E) ) 



(29) 



where (-) s = / °°(-)p s rf/?. On the other hand, since {(3) s — 
-(d\nB(E)/dE), 



dE 



d_ 
~dE 



B{E) jo 

W 2 )s + W)l 



pf((3)e-P E dp 



therefore we obtain a thermodynamical relation 
d 2 \nB{E) d(/3) s 



dE 2 



dE 



(30) 



(31) 



From Eq. H30JI , we obtain a discrete expression for B 1 {E) 
as follows. 
Putting (/3 2 } s 



e 2 , 



It follows that 
2edE 



dE 



1 



d(fl)> 



Solving this leads to 

09). 



-2nE 



d(fl)s 



W)s + ■ 



(32) 



(33) 



(34) 



where c is an integral constant. The assumption that 
the system's temperature on average would approach to 
zero if its energy vanishes by a classical description (i.e., 
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(f3) s — > oo when _E — > 0) provides c = 1. Then the 
solution becomes 



((3) s = -ecoth{eE). 



(35) 



We have 



lnB(B) = J ecoth(eE) = ln[a sinh(eS)] (36) 
where a is a constant of integration. Therefore 



JT 1 ^) 



a(e £ 
2e 



-eE 



ItE , -4e£ 



.(1 + e-^ + e 



(37) 



where £„ = 2(n + |)-\/ ((3 2 ) S E is discrete levels for the 
inverse of the superstatistics with the lowest -J ((3 2 ) S E. 
This discrete quantity has a dimension of specific heat, 
so that we can view that the system can be characterized 
by a superposition of many different levels of hierarchies 
with different classes of specific heat. However note that 
they have a common energy E. 



F. An inequality for B(E) 



where we have put the average as 

(•}c = £(-) e " BCc /£cOE), B c (E)=J2 e ~ E( °- ( 4 °) 

Under the condition that the deviation is small i.e., A£ <C 
1, 



(e~ EA <) c = 1-E{AQ C 



-E{Ai) c 



(41) 



Thus, within the first order perturbation, B(E) can be 
expressed as 



B(E) ~ B c {E)e- E ^-^)c 



(42) 



Furthermore, since an inequality (e^) > e^[28|. we can 
obtain an inequality which corresponds to the Gibbs- 
Bogoliubov-Feynmann inequality 29] in the usual statis- 
tical mechanics 



B(E) > B c (E)e- E{c -^c 



(43) 



It is worth noting that when £ is simply equivalent to 
(3 the equality holds and B(E) provides the Boltzmann 
factor, which corresponds to the environment with no 
fluctuation in (3. 



We consider a situation where single (3 in a Gibbsian 
statistical factor is replaced by more generic one. The 
sum of these factors which can be regarded as a discrete 
expression for the superstatistics may be written as 



B{E) 



{/3} 



(38) 



where C({/?}) is a quantity which is uniquely determined 
in a system once a distribution of /3 's is given. Let 
us consider a bound which is satisfied by B{E) from a 
perturbation theory about the distribution of the inverse 
temperatures over the system. We express the difference 
in C({/3}) between a certain realization of distribu- 
tion over a system and its coarse-grained one as A£({/3}) 
(i.e., AC({/3}) = C({/3})-Cc({/3}))- By the coarse-grained 
inverse temperature field we mean that the system has a 
single homogeneous one over the system and it appears to 
be reached at an equilibrium states if we look at whole 
system from far away. If we look closer at the system 
we should see the fluctuating field. This illustration also 
relates to the number of division of the whole system 
considered. The role of /Vs in the exponential function 
is just the same as the spin variables {s,} in the Ising 
Hamiltonian system, where the partition function and 
the system's energy are quantities for discussing the the- 
ory instead of B(E) and £ respectively. Therefore we can 
write 

B(E) = ^ e -E(UU3})+AC) = B c (E)(e- EA ^) c (39) 

{/3} 



III. DISCUSSION AND SUMMARY 

In this paper, we have investigated some thermody- 
namical structures and information aspects of the super- 
position of statistical factors which is achieved by the as- 
sumption of the Gibbsian local equilibrium distribution 
and the existence of the spacial fluctuation of the inverse 
temperature over the system. We spotlighted some struc- 
tural side of the superstatistics rather than the concrete 
application of it to the physical system. 

To emphasize the noticeable similarity in the structure 
to the canonical partition function, we can take a stand- 
point of the generating function[19|, which is constructed 
by a product of a statistical factor and a structure func- 
tion uu(x) (density of states). The generating function 
G{n) = / °° iv(n)e~ flu dfJ l can be related to the variance of 
its variables as 



G( M ) 



1 



1 



(44) 



where a 2 , = (v 2 ) — (v) 2 is a variance calculated from 
a definition (•) = J Q {-)uj{y)dv. \i and v represent ei- 
ther E or (3 respectively. When the symbols denote 
(//, v) — (E,(3) and (J3,E), the generating function rep- 
resents Z c {(3) and B(E) respectively. From this observa- 
tion, we arrived at a view that the canonical partition 
function and superstatistics have a universal relation to 
the degree of dispersion of environmental variables irre- 
spective of the form of density of states (for E and (3), 
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TABLE I: Gibbsian superstatistics and partition function 



Statistical factor 


Density of states 


Generating function 


Gibbsian 


P(E) 


Zc(P) 


e -PB 


m 


B{E) 



indicating that they are two different aspects of the sin- 
gle generating function (Tabled). In order to extract 
underlying structure of notion of superposition of statis- 
tical factors, we here study more generic circumstance 
for a statistical factor g(/3, E), which is assumed to be lo- 
cally realized. We write g(/3,E) = g((j3),E)^((/3},j3,E), 
where 0,E) = g- x ({f3),E)g{P,E). Therefore the 

superstatistics of the general statistical factor g becomes 



B(E)= / w(P)g(p,E)d0 = g((0),E){Z({0},0,E))- 
Jo 

(45) 

By Taylor expanding £((/3), /?, E) around (J3), we obtain 



B(E)=g(((3),E) 





H 


1+ r w 


/3=(/3> 



(46) 



This allows us to conclude that it is possible to relate pa- 
rameters contained in the distribution u>(/3) to the tem- 
perature dispersion a 2 for any statistical factors. Beck 
and Cohen showed that when u>{(3) is the gamma distri- 
bution and g = e~^ E an index q appearing in the Tsallis 
factor can be expressed as q — (j3 2 )/(/?) 2 @]. However, as 
we demonstrated above, this feature applies universally 
for the other combination of superposition of statistical 
factors |3l| once we assume the local equilibrium states 
for each partitioned region. We give two examples when 
the local equilibrium does not obey the Gibbsian form. 
In the cases of a Tsallis type power-low and a stretched 
exponential type local equilibrium, we respectively have 



B(E) 



-((S)E 



1 



1 



-qa 2 E' 



where e x q is the ^-exponential function frequently used in 
the nonextensive statistical mechanics y^] when crE is 
not large and 



q = 1 and rj = 1 because the statistical factors recover 
the Gibbsian. 

Although the superstatistics constructed in the sense of 
plays a role of a new distribution incorporated by sys- 
tem's global information on fluctuating inverse tempera- 
ture, it is not normalizable with respect to energy E. The 
situation is true for a partition function Z(f3) in the usual 
statistical mechanics, which gives / °° Z{[3)df3 ^ 1. We 
may construct a conjugate distribution function Eg. (|28fl . 
allowing us to write the field f(/3) = B(E)e" E p s for the 
purpose of explaining (fitting) observed distribution func- 
tions. A quantity In B~ 1 (E)/E can capture states of the 
system and can be finite even when E — ► oo|34| due to 
concavity of lnB~ 1 (E). In this paper, we did not con- 
sider the dynamical mechanism to keep the fluctuation 
/(/?) as a specific form. Instead, we would like to mention 
what we can say more about a relation between variance 
a and a thermodynamic force (affinity) when we regard 
Fit) = — In B(E)/E as a fundamental quantity to see 
the approaching process to uniform temperature over the 
system: 



AF(t) = F(t) - F 



E \B (E) 



where Bt(E) is a superstatistics calculated from f(j3) at 
time t and Bq(E) denotes one that from uniform temper- 
ature. Bo(E) is equivalent to the factor e~^ E . Then, 
in consideration of the universal correction term derived 
for all types of /(/?), up to the second order in energy, 
AF(t) is found to be AF(t) = -ln(l + \a 2 {t)E 2 ) /E, 

or equivalently a(t) — y/2(e~ EAF( - t ' ) — l)/E, expressing 
that the system has a single temperature in the course of 
the thermodynamics force vanishes. 

In summary, we have explored some formal proper- 
ties of the superposition of statistical factors proposed 
by Beck and Cohen. To understand what the notion 
of superstatistics represents, we have attempted to see 
links to information and thermodynamic relation. Non- 
equilibrium states or relaxation processes to uniform tem- 
perature field over the system would be interesting for the 
future research. 



B(E) 



where 2/(77) = -■qE-q(l3)^- 2 {-q - riE r >(/3) - 1). Then the 
universal quadratic correction terms to the alternative 
factors 1 + qa 2 E 2 /2 and 1 + y(ij)cr 2 /2 are respectively 
needed no matter what the precise form of the fluctua- 
tion is. Note that we have the same form as Eg . (|44|l when 
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